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We investigate the semileptonic decays of B and D mesons into vr and p mesons, 
respectively, by means of QCD sum rules. We find that for the vector form- 
factors involved the pole dominance hypothesis is valid to good accuracy with 
pole masses in the expected range. Pole dominance, however, does not apply 
to the axial formfactors which results in specific predictions for the predominant 
polarization of the p meson and the shape of the lepton spectrum. For the to- 
tal decay rates we find T{B^ ir+e~i>) = (5.1 ± 1.1) iKbP 10^^ s-\ T{D^ 
TT-e+u) = (8.0 ± 1.7)|FcdplO^°s-\ T{B° p+e-u) = (1.2 ± 0.4 ) lO^^ s'^ 

and r(L»o ^ p-e+i/) = (2.4 ± 0.7) 10^ s~K 
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I. INTRODUCTION 



Semileptonic weak decays of heavy mesons have proved to be a very important tool in 
exploring the Higgs-sector of the standard model and, in particular, the strength of weak 
decays of quarks, parametrized by the CKM matrix. Once the relevant hadronic matrix 
elements are known, the CKM matrix can be extracted from experimental measurements of 
the decay rates. Whereas the situation (both experimental and theoretical) now seems rather 
settled for the dominant decays B — > D^*^eiy and D Keu, at least at a level of accuracy 
of ~ 10% for the CKM matrix elements the experimental results for the Cabibbo- 

suppressed decays B, D ^ tc, peu at present still suffer from large statistical uncertainties. 
This situation, however, will improve in the near future, since the exploration of these decays 
is motivated by the quest for \Vub\ which of all the CKM matrix elements still is the one 
most poorly known 

There exist several theoretical calculations employing relativistic [^,[^ or non-relativistic 
P) quark models as well as (for D decays) lattice calculations P,p!0|, and, quite recently, 
attempts to relate formfactors of D decays to those of the corresponding B decays by means 
of the heavy quark effective theory [|lT|. All these models are, however, for conceptional or. 



as for lattice calculations, for economical reasons not capable of calculating the full dynamics 
of the decay process, even apart from model limitations. A quite standard procedure is to 
determine a formfactor / at some fixed point of t, the momentum transfer squared to the 
leptons, and then to assume either some pole-like t-dependence, 

fit) - (1.1) 

where rUpoi is the mass of the lowest lying resonance coupling to the corresponding current 
( 0, e.g.), or an exponentially increasing formfactor as in the non-relativistic model of 0. 
Indeed, at the level of a desired accuracy of, say, ~ 20% for the rates, the details of the 
functional dependence do not matter as long as tmax, the maximum value of t allowed by 
kinematics, is much smaller than rrip^i and the formfactors vary only slowly (as for B — >■ 
0^*^61^). A certain deviation from that insensitivity is noticeable in the decay D K*ey 
where all "conventional" (quark model and lattice) calculations are not capable to reproduce 
neither the absolute value of the rate nor the small value of the ratio of rates of longitudinal 
to transversal polarized K* (cf. |[I[]). The actual functional form of the t-dependence becomes 
crucial if t™'"'^ /m'^oi ~ 1 ^"^^ decays B ^ tt, pev which in the future will provide us 

with the most accurate information on Vub- Thus a point seems to be reached where an 
increased attention should be paid to the investigation of the t-dependence of formfactors. 

In fact, there is another method for calculating hadronic matrix-elements including non- 
perturbative effects which relies on the field-theoretical aspects and features of QCD and 
was designed to make maximum use of known manifestations of non-perturbative QCD, the 
QCD sum rules method Originally invented for the calculation of vacuum-to- meson 

transition amplitudes, it soon found application to the calculation of the electromagnetic 
formfactor of the pion |]13[ and other meson-to-meson transition amplitudes (cf. for 



a review). Although this method in general yields less detailed results than fine-tuned 
models, it has got the advantage that only a small number of parameters is needed that have 
an evident physical meaning (e.g. quark masses) and/or characterize the non-perturbative 



2 



regime of QCD (e.g. the so-called quark condensate, the order parameter of chiral symmetry 
breaking). Once these parameters are fixed from well known processes, they can be used to 
calculate for instance heavy meson decays. 

In previous publications |[l5| , p!6|JT7| , p!8[| , we have shown that the t-dependence of the form- 



factors oi D ^ and B D^*^ can reliably be calculated by means of QCD sum 

rules. As a general pattern to be followed, formfactors determined by vector currents ob- 
serve a pole-type behaviour with pole masses in the expected range, those determined by 
axialvector currents in general do not. In view of the above stated advantages of the QCD 
sum rules method, we feel it worthwile to apply it likewise to the Cabibbo-suppressed decays 
B, D ^ n, pev and to go beyond the existing QCD sum rule calculations p9|j20|j2l|j2^j23|j2^ 
which for reasons to be explained in the next section were restricted to a determination of 
the formfactors at t = 0. We expect to gain a reliable picture of the dynamics of these 
decays and well founded predictions of their decay rates within the scope of accuracy to be 
obtained by QCD sum rules, i.e. at the level of 20% at best. 

Our paper is organized as follows: in Sec. II we present the QCD sum rules method, 
improve the existing calculations and collect all necessary kinematics. In Sec. Ill we evaluate 
the sum rules and give results for the decays D ^ tt, pev and B ^ it, pev. Finally, in 
Sec. IV we discuss the results and compare it both to experiment and to other calculations. 
Formulae and technical details are collected in the appendices. 

II. THE METHOD 

We consider the three-point functions 

W = e j d'^xd'^ye-'P"'^+'^P"-P^'>y {0\Tjt''{0)Vi;i{y)jlj{x)\0) 

= t{PH+PLrpi^+ + --- (2.1) 

and 

T'^" = 1^ j d'^xd'^ye-'P^^^'^P^'P^^y {0\TjUO)iVhi - A;,;)'^(y)jl(a;)|0) 

= tg^'^To - t{pH + PLrPH^+ - ^^%yHPT^v + • • • , (2.2) 

respectively. Here {Vm — A^i)^ = —l^)h is the weak current mediating the weak decay 
of the heavy quark h with mass irih into the light quark / with mass m^, ju = Qi^bh is an 
interpolating field describing the pseudoscalar meson H {B or D) built up from h and the 
light antiquark q, and = 57*^(75)/ interpolates the light vector (pseudoscalar) meson 
L {it or p). ph and pi are the momenta of the heavy and the light meson, respectively. 
In the above equations, we have made explicit only those Lorentz-structures that actually 
contribute to the decays under consideration. 

The correlation functions (|2.1|) and ( p.2|) are functions of the scalars pj^, and t = 
{Ph — PlY can be calculated in perturbation theory for Euclidean values p\j — ml, 
p\ — mf <^ 0. On the other hand, the singularity structure of the correlation functions is 
known, and thus we can represent them by double dispersion relations in p'jj and e.g. 
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n+ = / dsndsL 7 — — — 1 w ^' ^ o \ + subtractions. (2.3) 
J [sh -Ph)[sl-Pl) 

The spectral function p^^^ can be expressed in terms of physical observables as 



rn.n 



m,n " i,j 



S {QL - EPLi)S {Qh - EPHj) 



{0\jnm){m\VM{n\jl,\0) (2.4) 

where one has to take the appropriate Lorentz-structure on the right hand side and qj^ = sh 
and g| = s^. The sum runs over all m- and n-particle states coupling to the currents jf^^ 
and jh-i respectively. In particular, we shall single out the ground states and write 

pfy^ {0\jt''\L){L\V/:,\H){H\jU0) + (2.5) 

where p'^""'^ contains both the contributions of higher resonances with appropriate quantum 
numbers and of many-particle states. The first term on the right hand side contains exactly 
the quantities we are interested in: 

{n\V>:i\H) = U{t){PH+P.), + f-{t){pH-p.),, (2.6) 
(P, mi - Ki\H) = -t{mH + m,)/li(t)ef ) + '^f^ {e*^'^PH){pH + P,) 

+ ^^^{e<%H){PH - P,) + (2.7) 

These are the relevant matrix-elements governing the hadronic part of the decays in question, 
decomposed in terms of the formfactors f±, and V, where uih and nip are the masses 
of the H and the p meson, respectively; A denotes the polarization state of the p. In the 
hmit of vanishing lepton mass, the formfactors /_ and do not contribute to the decay 
rates and henceforth will not be considered. Expressed in terms of the above formfactors, 
the spectra with respect to the electron energy E read: 



dr{H n+e-p) G 



dE 167r^mij 





■"max 

j dt{2E{m]j -ml + t)- m„{t + 4E^)}fl{t), (2.8) 



dr{H p+e-p) 



dE 



1287r3mi 

^ 



'■max 

j dtt{{l-cosef H^_ + {1 + cos ef Hi + 2{l-cos^e) H^} (2.9) 



with the helicity amplitudes 
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AV2 

H± = {mn + mp)Ai{t) t ■ ^W, (2.10) 

mn + mp 

1 f A^/^ 1 

^0 = 1= « -ml- t){mH + mp)A,{t) ^(t) . (2.11) 

2mpVt [ ^ niH + mp J 

tmax, the maximum value of t, the invariant mass squared of the lepton-pair, is given by 

niH ^ . (2.12) 

niH — It, J 

9 is the angle between the p and the charged lepton in the (e"z/) CM system and given by 

cos ^ = ^ {m]j -ml + t- AttihE) (2.13) 
V A 

where A = {mjj + — t)^ — Am^ml. 

Returning to ( p.l|) and (|2.2|) , it was the idea of Shifman, Vainshtein and Zakharov [|1^ to 



account for non-perturbative corrections to correlation functions by expressing them via an 
operator product expansion (OPE) including terms that vanish in the perturbative vacuum, 
but acquire finite values in the QCD vacuum. These so-called condensates characterize the 
long-distance behaviour of the correlation function and we are led to write, e.g. 

n+(p^,pi, ^) = E ^i\p'H.pl t) ( 1 a 1 ). (2.14) 

n 

The Wilson coefficients 11^'' can be calculated with the aid of perturbation theory for nega- 
tive values of pjj^ — ml and p| — mf. The ( | 0„ | ) are vacuum expectation values of gauge 
invariant operators, the above-mentioned condensates (cf. App. A). The first term in the 
series just covers usual perturbation theory, the others are non-perturbative corrections. In 
our analysis we will take into account the lowest dimensional condensates up to dimension 
6, where we improve existing calculations [25|| by the inclusion of the contribution of the 



gluon condensate (App. A). Equating ( |2.3| ) and the OPE ( |2.14| ) yields expressions for the 
formfactors determining ( p.l|) and ( p.2|) in terms of QCD parameters (like quark masses) and 
condensates. Before, however, we can start to evaluate these sum rules, we have to specify 
how to treat p'^°'^^ in (|2.5|) . As for that, we employ the argument of quark-hadron-duality 
P^ , e.g.) and model by the contribution of usual perturbation theory above some 



thresholds and s° , 



^cont ^ ^pcrt _ 0^^O _ 0(^0^ _ (2.15) 

which restricts the region of integration in the plane of sh and sl to the gray area depicted 
in Fig. m. The calculation of p^'^'^^ for t > involves some delicate points connected with the 
possibility of the spectral function to become singular. For the discussion of the additional 



"non-Landau" contributions caused by these singularities we refer to [T^. In the numerical 
analysis we will tacitly include those contributions whenever necessary. 

The dependence of the sum rule on the continuum model as well as the error induced 
by truncating the OPE series can be diminished by the application of a Borel transfor- 
mation. For an arbitrary function of Euclidean momentum, /(-P^) with = — p^, that 
transformation is defined by 
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1 

f:=BMM')f= lim TTji'PY^' 77^^^ (2-16) 

P^/N=M^ fixed 



where is a new variable, called Borel parameter. For a typical term appearing in the 
OPE, the transformation yields 

^-(^'^) = j^. (-1)" p^i.^-"'""- (2-1^) 

Since condensates with high dimension get multiplied by high powers of (p^ — m^) in the 
denominator, their contributions get suppressed by factorials. In addition, the contribution 
of higher resonances and the continuum, p™"^*^ gets exponentially suppressed relatively to 
the contribution of the ground state, which is just the desired effect. 

We are now in a position to write down sum rules for the relevant formfactors. From 



3) and we find: 



Parts of the explicit formula for the correlation functions can be found in [jl5|. For the 
present analysis, we in addition have calculated the contributions of the gluon condensate 
and the contribution of the four-quark condensate to /+; the formulae can be found in the 
appendices. Note that we have expressed the vacuum-to-meson transition amplitudes in 
terms of the corresponding leptonic decay constants as 



(0|t/^75&|i?°) = /B — , (2.22) 
nib 

{Q\d^,u\p-^,\) = fpm/^\ (2.23) 

( 0| J7^75M 1 7r+ ) = iUp^y. (2.24) 
The question of how to treat these quantities will be discussed in the next section. 



III. EVALUATION OF THE SUM RULES 



In the numerical evaluation of the sum rules ( |2.18|) to ( |2.21|) we use the following values 
of the condensates at a renormalization scale of 1 GeV: 
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(gg)(lGeV) = (-0.24GeV)^ 



(^G') = 0.012 GeV^ 



TT 



{qagGq){l GeV) = 0.8 GeV^(gg)(l GeV), 



16 

u,d,s 
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4nas{duud) ^ Anasiqq)'^. (3-1) 

We use leading-order anomalous dimensions of the quark and the mixed condensate to 
evaluate them at a scale /i which is given by the harmonic mean of the Borel parameters, 



fi^ = yM^M^. For the four-quark condensates we assume vacuum saturation. Since their 
contributions are tiny, we neglect the scale-dependence. In general the sensitivity of the 
formfactors on the actual values of the condensates will be smaller than 10% when changing 
{—{qq)Y^^ by 10 MeV and most pronounced for the axial formfactors A2. The smallness of 
the contributions of the contributions of the four-quark condensates indicates that higher 
order power corrections are well under control. Concerning quark masses, we put the masses 
of the u and the d quark to zero, for the heavy quarks we use the renormalization-group and 
-scheme invariant pole mass. Its connection to the running mass in the MS scheme is given 
by (for scales fi <^ '^Ws) 

mpoie = m^{^^) + ^ + 1^ ^ ] ■ (3-2) 



The numerical values are (cf. [|14|; a recent determination of is given in |p6|| ) 

mfo = (4.6-4.8) GeV, = (1.3-1.4) GeV. (3.3) 

For the leptonic decay constants we use the experimental values = 0.133 GeV and fp = 
0.216 GeV. For fs and fo we employ two-point sum rules ( e.g.), discarding radiative 
corrections. We expect the accuracy of the sum rules for the formfactors to be increased by 
that, since both in the limit of infinitely heavy quarks and for the matrix-element {B\Vfj,\B), 
where charge conservation fixes the formfactor at zero recoil, QCD sum rules yield the 
correct normalization independent of the values of quark masses, continuum thresholds and 
Borel parameter ||T7| , ^ , ^ , provided the continuum thresholds in the two-point and the 
three-point sum rules are chosen equal and the Borel parameter in the three-point sum rule 
takes twice the value of that of the two-point sum rule. We will take these prescriptions over 
to the case where the outgoing meson is light, and actually the sensitivity of the resulting 
sum rules on rrib or rric is greatly reduced as compared to the sum rule for fs or fjj. In 
addition, the effect of the unknown radiative corrections to the three-point function should 
tend to cancel against the radiative corrections to the two-point function. Concluding, we 
take both the range of Borel parameters, the "sum rule window", and the values of the 
continuum thresholds from the two-point sum rules, i.e. we evaluate ( |2.18| ) to ( |2.21| ) in the 



range 7 GeV^ < < 10 GeV^ and 2 GeV^ < < 4 GeV^ and for continuum thresholds 
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s% = (6-7) GeV^ 4 = (34-36) GeV^ s° = (0.75 - 1.0) GeV^ s° = (1.25 - 1.5) GeV^ In 
addition, we choose a fixed ratio of tlie Borel parameters, 

Tliis procedure ensures tliat perturbative and non-perturbative corrections in botfi tfie fieavy 
and tlie liglit cfiannel are equally weighted. The sum rules are rather insensitive to the actual 
value of that ratio, and changing for example M^/M^ from 3 to 5 results in changes of the 
results of at most 10%. 

Our aim is to extract the t-dependence of the formfactors from the sum rules, and 
thus we are restricted to a range of values where the correlation function can be expected 
to be reliable in that variable. That is, we have to stay approximately 1 GeV^ below 
the perturbative cut starting at t = rn^^. Thus we can trust the sum rules up to t ~ 
20 GeV^ for the 5-decays and t ^ 0.9 GeV^ for the D-decays. The maximum values allowed 
by kinematics are t^^/ = 26.4 GeV^ = 20.3 GeV^ t^^/ = 3.0 GeV^ and tgZ' = 

1.2 GeV^, so apart from i? — > p we cannot cover the full range of t. Albeit the accessible 
range is sufficient to determine the shape of the formfactor and the total rate, it is not 
for the calculation of the electron spectrum. We thus take the attitude to extrapolate the 
formfactors to tmax to get the electron spectrum for large values of the electron energy E. 
This procedure does not introduce too large an uncertainty for D ^ p since according to 
the above remarks at least 80% of the integration range in t is covered, and for D ^ tt 
and B ^ TT we will find a pole-type behaviour of the formfactors which facilitates the 
applicability of the extrapolation. 

In Fig. 0(a) the formfactor f^^^{t = 0) is shown as function of the Borel parameter 
M^. The different curves correspond to different choices of the set of input parameters. To 
be specific, we use = 1.3 GeV, s% = 6 GeV^ s° = 0.75 GeV^ (set CI), = 1.3 GeV, 
s% = 6GeV^ 5° = IGeV^ (set C2), = 1.4 GeV, s% = 7GeV^ s° = 0.75 GeV^ (set 
C3), m, = 1.4 GeV, s% = 7GeV^ s° = 1 GeV^ (set C4). The value of s% is taken as the 
best-fit continuum threshold for the sum rules for f^, s° is taken from [|12|. In the "sum 



rule window" 2 GeV^ < < 4GeV^ the f^^'^{0) is quite stable and the dependence on 
the values of both the mass of the c-quark and the continuum thresholds as respresented by 
the spread of curves is well under control, yielding /+(0) = 0.5 ±0.1 where the error is an 
educated guess based on both the dependence of the sum rule on the input parameters and 
the intrinsic uncertainty of the whole method. Perturbation theory and quark condensate 
give the dominant contribution, the other condensates contributing at the level of ~10%. 
Thus the series of power corrections is well under control. 

/+^'^(0) as function of the Borel parameter is shown in Fig. 0(b). Here we use the 
parameter sets Bl (mf, = 4.6 GeV, s% = 36 GeV^ s° = 0.75 GeV^), B2 {rrib = 4.8 GeV, 
4 = 36GeV^ si = IGeV^), B3 {nib = 4.8 GeV, s% = 34GeV^ s° = 0.75 GeV^), B4 
(mfe = 4.8 GeV, s% = 34GeV^, s° = 1 GeV^). Again the value is remarkably stable against 
variation in the quark mass, in the continuum thresholds and the Borel parameter. Form 
Fig. 0(b) we find /+^'^(0) = 0.26 ±0.02. This value is higher than obtained in which is 
due to the contribution of the gluon condensate not included there. 

In Fig. ^(a) we show f^^^{t) as function of t, normalized to its value at t = 0, which 
representation emphasizes the differences in shape. We have chosen = 3 GeV^ and give 
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curves for all parameter sets. We find a rise in t which is very well compatible with a pole- 
type behaviour as suggested by the pole dominance hypothesis ( |1.1|) . From a pole fit we get 
fT^poi = (1-95 ± 0.10) GeV (including all sets and Borel parameters within the window). If 
pole-dominance were exactly valid, the pole-mass would be rriD* = 2.01 GeV, so QCD sum 
rules confirm pole-dominance for D ^ n. 

Pole-dominance is likewise valid for the B —>■ tc transition whose normalized formfactor 
is depicted in Fig. ^(b) as function of t and for all parameter sets at = 8GeV^. Pole 
fits yield pole-masses of ~5.1 GeV for the sets Bl and B2 and ~5.2 GeV for B3 and B4. 
From that we are forced to exclude the lower value of the b-quark mass, m^, = 4.6 GeV, 
from our analysis (since the formfactor would become singular at tmax) and stick to rrib = 
4.8 GeV. The "physical" pole is at rriB* = 5.33 GeV what nicely agrees with the fit value 
(5.25 ± 0.10) GeV from sets B3 and B4. 

In Fig. I we show the electron spectra dT/dE as functions of the electron energy E which 
can be obtained from the formfactors f^^'^ and /+^'^, extrapolated up to tmax according 
to pole- dominance. 

Let us now turn turn to the decays H —>■ pev. In Fig. ^ we show the formfactors of D —>■ p 
at t = as functions of the Borel parameter for the parameter sets C5 {rric = 1.3 GeV, 
s% = 6 GeV^ s° = 1.25 GeV^), C6 (m, = 1.3 GeV, s% = Q GeV^ s° = 1.5 GeV^), C7 {m^ = 
1.4 GeV, s% = 7GeV^ s° = 1.25 GeV^), C8 (m, = 1.4 GeV, s% = 7GeV^ = 1.5 GeV^). 
All formfactors are quite stable and we find A^^p{0) = 0.5 ± 0.2, ^^(0) = 0.4 ± 0.1 and 
V^^p[0) = 1.0 ± 0.2 where as in the previous cases the error is intended to include likewise 
systematic uncertainties. 

For 5 ^ p we find from Fig. | Af^''(0) = 0.5 ± 0.1, A^^^iO) = 0.4 ± 0.2, ^^^''(0) = 
0.6 ± 0.2 with the parameter sets B5 (m^ = 4.6 GeV, s% = 36 GeV^ s° = 1.25 GeV^), 
B6 {nib = 4.6 GeV, s% = 36 GeV^ s° = 1.5 GeV^), B7 {nib = 4.8 GeV, 4 = 34GeV^ 
s° = 1.25 GeV^), B8 {nib = 4.8 GeV, s% = 34GeV^ s° = 1.5 GeV^). All these formfactors 
depend only slightly on quark masses and continuum thresholds and are stable in the Borel 
parameter, except for ^^^''(0). Here we observe for B7 and B8 a rather strong dependence 
on the reason being the extremely small contribution of perturbation theory which is of 
about 10% only. 

In Fig. 1^ the normalized formfactors of the D p transition are plotted as functions of 
t for = 3 GeV^ and all sets of parameters. We find a decrease of A^^^{t) in t which is 
nearly independent of the parameter set used. This behaviour is in clear contradiction with 
pole-dominance which predicts an increase determined by the pole-mass m^jn- = 2.42 GeV 
corresponding to A^^^{t)/A^^^{0) = 1.21. A similar behaviour is encountered for A^^'' 
where we find a decrease of about 10% at t = 1 GeV^ depending on the parameter set 
used. For the vector formfactor we have an increase in t with a best-fit pole-mass of nipoi = 
(2.5 ± 0.2) GeV which is a little bit larger than predicted by pole- dominance. 

For the normalized formfactor Af^''{t)/Ai^''{0), shown in Fig. ^(a) at = 8 GeV^ 
for all parameter sets, we find a rather unexpected shape with a minimum at t ~ 15 GeV^. 
Formally, this minimum is due to the interplay between decreasing contributions of pertur- 
bation theory and quark condensate and an increasing one of the gluon condensate which 
becomes effective at large t. For A^~*''(t)/A^^''(0) (Fig. |^(b)) we find a moderate increase in 
t which at large t is again compensated by a negative contribution of the gluon condensate. 
For V^-'P{t)/V^^f'{0) Fig. |(c) shows the usual increase in t corresponding to a pole mass 
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of (6.6 ± 0.6) GeV, about 1 GeV larger than predicted by pole-dominance. 

Finally, in Fig. ^ we show the electron-spectrum dT/dE of the decay D — > peu as function 
of the electron energy E for set C5 and = 3 GeV^ where the formfactors are extrapolated 
in the range t > 1 GeV^. Fig. |l^(a) shows the electron spectrum dT/dE of I? — peu as 
function of E for set B7 and = 8 GeV^ as quite sharp and concentrated around large 
electron energies. 



IV. RESULTS AND DISCUSSION 

In the previous section we have given a careful analysis of the semileptonic heavy-light 
decays D — >• vrez/, B —>■ neu, D —>■ peu and B — > peu. We have put some stress on 
the calculation of the t-dependence of the formfactors which for the vector formfactors in 
general can well be described by a pole- dominance formula, whereas the axial formfactors 
tend to decrease in t and even develop extrema. The numerical results of our calculation 
as well as other models are collected in the tables. The formfactors at t = can be found 
in Tables I and II, the rates in Tables III and IV. The rates were calculated either using 
pole-dominance (as indicated in the tablenotes) or some other model for the t-dependence. 
In addition to the total rates we give for H — > peu the ratios Tl/Tt and r+/r_ where 
the index denotes the polarization state of the p (longitudinal, transversal, positive, and 
negative helicity, respectively). The corresponding electron spectra are shown in Figs. § 
{D, B vrez/), g{D ^ peu) and |TO|(a) {B peu). 

The only decay, where a comparison to experiment is possible so far, is D —>■ neu. In 
addition there exist several model calculations in literature, using QCD sum rules 



quark models f^J^J^ and some lattice calculations |]9|JT0[|. One calculation relying on the 
heavy quark effective theory |11[] takes the experimental result as input to their values 
of the formfactors of 5 — > tt, peu. The theoretical predicitions of r{D vrez/) differ by a 
factor of two, and assuming \Vcd\ = 0.22, which can be inferred from the unitarity of the 
CKM-matrix with high accuracy, we find that the central value of our rate is by two standard 
deviations smaller than the experimental value. This discrepancy is not strong enough to 
be conclusive and might be due to the neglection of radiative corrections to our sum sules. 
Still further experimental effort in improving statistics is to be desired to clarify this point. 



For D pev experiment only has set an upper bound for the total rate so far |^ 
Our value is by a factor five smaller than these. For the ratio T{D p)/T{D tt) we 
get 0.3 which again is smaller than the predictions in other models which yield a maximum 
values of 1.8 §]. We remind that the corresponding ratio for the Cabibbo favoured decays, 
r(_D K*)/V{D K), is approximately 0.5 and that we do not expect fiavour SU(3) 
to be broken by a factor of two or more. The formfactors at t = roughly agree in all 
models except for the lattice calculation [|l^ which predicts vanishing A^^'^^O) and a small 
value of V^^~*''(0) yielding a large value of Tl/Tt- 

For the b ^ u decays we do not dare to quote any experimental upper bound for the total 
rates due to the uncertainty in \ Vub\ (but cf. 0). We remark that the total rates for B —>■ Treu 
summarized in Table IV and obtained by QCD sum rules [p2| , p3tp4| , the quark models 



and the HQET calculation [|TT| differ by a factor of 26. For B peu this value shrinks to 4. 
That spread in predictions clearly shows the neccesity for an accurate investigation of the 
t-dependence of the formfactors that we have concentrated on in this paper. With the t- 
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dependence obtained by QCD sum rules we obtain T{B^ p'^ev) = (1.2±0.4) |T46p-10^^s ^ 
where the p has mainly negative helicity. Furthermore, we find r{B p)/T{B ^ vr) = 2.4 



which is smaller than all other model predictions ranging from 3.1 to 11 except for |11 



which predicts 0.6. In Fig. p!U|(b) we give the electron spectrum 1/V dV / dE for B — >■ pez/ as 



obtained in this paper (same parameters as in Fig. |TO|(a)), in the BWS model (using pole- 
dominance) and the non-relativistic GISW model . The chosen normalization emphasizes 
the difference in shape rather than in the absolute normalization. Although the BWS 
spectrum is softer in the endpoint region above the threshold for charm-production, the 
only region where 6 — > m transitions can be observed, our spectrum and that of GISW are 
nearly indistinguishable for E < 2.4 GeV. If, however, a detection of the polarization of the 
p was feasible, one could test the considerable deviations of the corresponding spectra in 
the different models. We predict the p to have predominantly negative helicity (as indicated 
by the very small values of Tl/Tt and r+/r_) whereas in other models the ratio Tl/Tt is 
closer to one. 

In Fig. |ll| we give a comparison of the inclusive h ^ u semileptonic spectrum calculated 
by means of QCD sum rules in [^] to the exclusive decay spectra B neu and B 



peu calculated with the same parameters {mb = 4.8 GeV, s% = 34GeV^ = SGeV^). 
Fig. p!T|(a) shows the spectrum in the restframe of the decaying B-meson, Fig. 0(b) in the 
laboratory system of an e~^e~ collider operating on the T{4S) resonance. From both we find 
that at high electron energies B pev constitutes nearly the whole differential inclusive 
rate, so it is worthwile to concentrate on measurements of the exclusive channels, where 
theoretical predictions are still not at their best precision, but are much more accurate than 
calculations of the inclusive spectrum (cf. 
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APPENDIX A: THE WILSON COEFFICIENT OF THE GLUON CONDENSATE 

In the following we present a technique for calculating Borel transformed Wilson coef- 
ficients directly from the loop-integrals. This method does not allow for the subtraction 
of continuum contributions, which, however, does no harm in our case as the total contri- 
bution of the gluon condensate to the three-point sum rule is small by itself (^ C(10%)), 
and so is its continuum portion. Besides, one would expect typical continuum contributions 
to show up as incomplete Gamma functions in the Wilson coefficiens, which, however, are 
absent in our formulae (i.e. in the sum of all diagrams, but are encountered in each diagram 
separately). Thus one is led to conclude that those contributions are actually absent in the 
processes under consideration. 

We calculate the diagrams shown in Fig. ^ that contribute to the Wilson-coefficient of 
the gluon condensate (asG^/vr) in the fixed point gauge 

x^'Al{x)=Q (Al) 
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with the gluon field A^, a G {1, 2, ... , 8}. For massless Ught quarks and with the coordinates 
chosen as indicated in the first diagram, diagrams I and II evaluate to zero. Note, that for 
massless light quarks there is no mixing of the gluon with the quark condensate. 

In the evaluation of the remaining diagrams we encounter integrals of type (since the 
Borel transform removes UV divergencies, there is no need for dimensional regularization of 
these divergencies and we thus stay with four-dimensional integrals) 



Although the sum of all diagrams is IR convergent, IR divergent terms occur at each step of 
the calculation and need proper regularization. We take the attitude to let the mass of the 
u-quark finite in the denominator of its propagator, > (but let m„ = in the traces) 
and regularize the singularities in the q-quark line (for a = 2) by shifting the power of the 
q-quark propagator to 2 — e (which from a technical point of view is simpler than introducing 
dimensional regularization). This procedure has got the advantage that all integrals with 
6, c > 1 can be obtained from the case 6 = c = 1 by taking derivatives with respect to the 
quark masses: 

^ c) = d{miy-^ d{mlY-' ^mM2-/.n(«, 1, 1)- (A3) 

Continuing to Euclidean space-time and employing the Schwinger representation for prop- 
agators, 

1 1 °° 

[p2 + ^2]a r(a) y ' ^ ' 

we find for the scalar integral n = with a — h — c — 1 (with capital letters denoting 
Euclidean momenta): 

/(1, 1, 1) = -ijdad(3dj exp(-Si^^ - - -^T- jml) (A5) 

where 

K^K + ^{(3P, + ^Pb), (A6a) 

E = a + /3 + 7, (A6b) 
T = -t. (A6c) 
The above representation proves very convenient for applying the Borel transformation with 

Bp2{M^) e""^' = 5(1 - aM^). (A7) 

From that, we get 
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7(1, 1,1):= Bp.{M'J Bp.{Ml)I{l, 1, 1) 



I 

m2 + M2 



1 ^-tmi+Mi) (A8) 



'6 

where the exponential integral function is given by 



oo 

Ei{x) = - Jdt ^ (A9) 



-X 



and 

ml t 

Actually IR divergent diagrams only occur for the scalar case where we find 



,2 ,,2 



7(2 - e, 1, 1) = — ^ ^ ' <; - + 1 - 27b - In -At - At 1 - Inz J> . (All) 

Here // is some arbitrary scale introduced to render the canonical dimension of the integral, 
which, however, cancels in the complete expressions for Wilson-coefficients, as it should. 
For larger values of n, we get 

t n 1 iv_ ( ^'^^Ml y 1 (p,,, PB,A 

(^^-mllMl-mllMl ^ ^-^/(M^^+M^) + ^Ei(-z))) , (Al2) 

7^,^, (a < 2, 1, 1) = Y(^ VM2 + M2 J Mp7| ~ 

+ ^1 W2t^(3-a,0;z) 



2(3 - a) \Ml Ml 



M^M^ 



C/(4-a,l;^)|, (A14) 

where we have continued back to Minkowski-space. U is the confluent hypergeometric 
function defined as 
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1 °° 



U{t,r,x) = — - Idtil + ty-'-'f-^e-^'. (A15) 
V{i) 



In addition, we use 



= [MT^MT' M-HM,r^W-..(«,^c). (A16) 

We now decompose the Lorentz-invariants A occuring in the Borel transformed correlation 
functions (|2.1|) and ( p.2| ) as 

A = J2k^^\On) (A17) 

n 

where (Oi) = (1) = 1, {Os) = {qq), {O,) = (a.^Vvr), {O,) = (qagGq), {O^,) = 
T!'as{qj'^X^qJ2u,d,sQlT^^q) and {Oq^) = 4:nas{duud) are the condensates taken into account. 
The formula for n G {1, 3, 5, 6} can be found in |T^ (with the same notations), where vac- 
uum saturation for the condensates with dimension 6 is assumed; the formulae for n = 4 are 
new and read for the relevant invariants: 



ii I - 



1 ^ 1 mlMUMi{Mi + 3M'J-ml{Mi + M'J} 



96M^ Ami 2AMl{M^ + Mlf'. 



9ml [M^ + Ml) 



2^2 



4M,^(M,^ + Ml) - ml{2M^, + Ml) mfM^ 



A8M^{M^ + Ml)z 16Ml{M^ + MD^z 



(A18) 



f(4) 



° M2M2 



m 



Ml , m^M„HM2(3M2 + lOM^) - Sm^M^} 



+ 



AM^ (Ml + Ml)^ IQM^{MI + ) 2 z^ 



mlMl{mlMi{23Ml + TMj) - 2M^{l?,Ml + SM^) - Sm^M^} 
MliSM^MliMi + 2MI) + AmlMliM^ - AM^Ml - UMl) - miMl} 



+ 



9QMl{Ml + M, 



2-^2^2 



mtMl{7Mi + 16Ml) Ml{AMi{2Ml - Mj) + ml{2Mi - 9Ml)} 
96Ml{Ml + Ml)^z'^ A8M^ {Ml + M^)z 

mlMt Mi{8Ml-5Mi)-3mlMl (Mj + Ml)z 

32M6(M2 + M^)z 96M^ 96M^ 



(A19) 



^(4) ^ mke-</^S 
MiMl 



1 ^ mlMliMi + ml) , M2(4M2 - ml 



96Mi 2AM^{Mi + MD^z^ Am^{M^ + Ml)z 



14 



+ 



16M6(M2 + M2) 



2^4 



(A20) 



,(4) _ rribC 



-ml/Ml 



mtMt 



V 



+ 



+ 



M2M2 [8M6(M2 + M2)224 12M6(M2 + M2)2z3 

M2{8M,4M2-4m2M2(M2 + 2M2)+m^M2} , 1 



48M6(M2 + M2)-z 



2^2 



48M2 



2M^{2M^-Mi 



24M^{Mi + M^Jz 



(A21) 



Note that we have checked our method for calculating Borel-transformed Wilson-coefficients 
for the matrix element {B\Vfj,\B) at t = where the result is uniquely determined by 
charge conservation. 



APPENDIX B: THE WILSON COEFFICIENT OF THE 
FOUR-QUARK CONDENSATE 

In addition to the formulae given in |T5[ , we also have calculated the contributions of the 



four-quark condensate to the decays B, D ^ neu which read (in the notation of (|A17|) : 



^(6,) _ m.e-^'l/^' ( 1 ml I t ml - 1 



M2M2 [9M^^9M^^{ml-t))- ^ ' 
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TABLES 

TABLE I. The formfactors of the c ^ d transitions at t = in different models. 



Reference /f^^ ^ ^2 V'^^p 

This paper 0.5±0.1 0.5±0.2 0.4±0.1 1.0±0.2 

I!" 0.7±0.2 _ _ _ 

@" 0.75±0.05 _ _ _ 

0.69 0.78 0.92 1.23 

0.51 0.59 0.23 1.34 

§^ 0.58±0.09 0.45±0.04 0.02±0.26 0.78±0.12 

0'=''^ 0.84±0.12±0.35 0.65±0.15±1[]:|| 0.59±0.31±lg;| 1.07±0.49±0.35 

0.79 0.55 0.28 1.01 
m 0.80l°:?l 



QCD Sum Rules 
^ Quark Model 
^ Lattice Calculation 

First error statistical, second systematical. 
^ HQET + chiral perturbation theory; value of /+ taken from experiment. 
f Experiment (using pole dominance) 

TABLE II. The formfactors of the h ^ u transitions at t = in different models. 



Reference /f^^ A'^^'' ^ F^^^ 

This paper 0.26±0.02 0.5±0.1 0.4±0.2 0.6±0.2 

H"''' 0.26±0.01 _ _ _ 

H"''' - 0.96±0.15 1.21±0.18 1.27±0.12 

0.23±0.02 0.35±0.16 0.42±0.12 0.47±0.14 

0.4±0.1 _ _ _ 

[g'^ 0.33 0.28 0.28 0.33 

0.09 0.05 0.02 0.27 

Ol'^ 0.89 0.21 0.20 1.04 



QCD Sum Rules 

^ Analysis suffering from a missing factor 12 in the perturbative contribution. 
^ Quark Model 
HQET + chiral perturbation theory. 
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TABLE III. Decay rates of the c ^ d transitions in units | Vc^p 10^^ s~^. denotes the portion 
of the rate with a longitudinal polarized p, with a transversely polarized p, r_|_ with a p with 
positive, r_ with a p with negative helicity. 



Reference 


r(i:)0 ^ TT-e+iy) 


T{D^ p-e+v) 


Tl/Tt 


r+/r_ 


This paper 


0.80±0.17 


0.24±0.07 


1.31±0.11 


0.24±0.03 





■■-•^^-o.ri 


- 


- 


- 




±.DD_Q 21 






- 






1 /i4_i n 


n Q 
u.y 






1.41 


1.38 


0.91 


0.19 




1.41 


1.40 


0.80 


0.13 


i 


0.77 


1.35 


1.33 


0.11 






0.83±0.19 


1.86±0.56 


0.16 




2m1:lit 


1.09 


1.10 


0.18 




1.9 


0.93 


1.40 


0.14 




|0.22/y,rf|2 • [i.dtli) 


l/\V,d?T{D+^ 


p°e+i/) < 0.71 • 10" s" 


-1 



" Rates calculated using pole dominance with mi- 

^ No values of formf actors given. 
Values of formfactors at t = identical to 
Values without errors from central values of Table I. 

^ Value of T{D^ — > 7r~e^z/) taken from experiment. 

f Experiment 



2.01 GeV, mi+ = 2.42 GeV. 



TABLE IV. Decay rates of the b ^ u transitions in units iKibP 10^'^ s~^. denotes the portion 
of the rate with a longitudinal polarized p, Tt with a transversely polarized p, with a p with 
positive, r_ with a p with negative helicity. 



Reference T{B^ n+e-u) TjB^ p+e-j?) Tl/Tt r+/r„ 

This paper 0.51±0.11 1.2±0.4 0.06±0.02 0.007±0.004 

0.68±0.23 - - - 

H - 0.77±0.42 

11" 0.302± 0.005 3.3±0.3 0.88l[j:^|^ 0.121°;°^ 

11^ 1.45±0.59 - - - 

[§° 0.74 2.6 1.34 0.16 

0.74 2.30 0.54 0.02 

§ 0.21 1.63 0.75 0.08 

0'^ 5.4 3.4 0.36 0.14 



Rates calculated using pole dominance with m^- = 5.33 GeV, mi+ = 5.71 GeV. 
Rate calculated using a modified pole dominance with m^- = 5.33 GeV. 
Values of formfactors at t = identical to . 
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FIGURES 



FIG. 1. The integration region for the perturbative contributions to the sum rules in the sl-sh 
plane (without non-Landau contributions). 

FIG. 2. (a) The formfactor /f^~"^(0) as function of the Borel parameter M^. The parameter 
sets in the legend are rric = 1.3 GeV, s% = 6GeV^ s° = 0.75 GeV^ (set CI), m,,, = 1.3 GcV, 
s% = 6GeV2, 4 = IGeV^ (set C2), rUc = 1.4 GeV, s% = TGeV^, 4 = 0.75 GeV^ (set C3), 
mc = 1.4GeV, s% = TGeV^, 4 = IGeV^ (set C4). (b) /f^''(0) as function of the Borel 
parameter with the parameter sets Bl (m^ = 4.6 GeV, s% = 36GeV^, s° = 0.75 GeV^), B2 
{nib = 4.8 GeV, 4 = 36GeV2, 4 = IGeV^), B3 (m^ = 4.8 GeV, 4 = 34GeV2, 4 = 0.75 GeV^), 
B4 {rub = 4.8 GeV, 4 = 34GeV2, 4 = 1 GeV^). 



FIG. 3. (a) The formfactor f^^'^{t)/f^^'^{0) as function of t for all parameter sets at 
M2 = 3GeV2. A pole fit yields n^pol = (1.95 ± 0.10) GeV, the maximum physical value of t is 
^max ^ 2.96 GeV^. (b) Like (a) for /f ^'^(t)//f "^'^(0). nipoi = (5.25 ± 0.10) GeV (for B3 and B4), 
^max ^ 26.4 GeV^ 

FIG. 4. The electron spectra dV/dE as function of the electron energy E for (a) the decay 
D irev (set CI, = SGeV^), (b) B irev (set B3, = 8GeV2). 

FIG. 5. The formfactors of D ^ p at t = as functions of the Borelparamctcr M^. The param- 
eter sets are C5 (m^ = 1.3 GeV, 4 = GGeV^, ^ = 1.25 GcV^), C6 {m^ = 1.3 GeV, 4 = 6GeV2, 
4 = 1.5 GeV^), C7 {nic = 1.4 GcV, 4 = 7GcV2, 4 = 1.25 GeV^), C8 (mc = 1.4 GeV, 
4 = 7GeV2, 4 = 1.5GeV2). (a) Af^^(O), (b) ^^^^^(O), (c) ^^^^(0). 

FIG. 6. The formfactors of -B ^ p at i = as functions of the Borelparameter M^. The param- 
eter sets are B5 {nib = 4.6 GeV, s% = 36GeV^ 4 = 1.25 GeV^), B6 (m^ = 4.6 GeV, 4 = SeCeV^, 
4 = 1.5 GeV^), B7 (m^ = 4.8 GeV, 4 = 34GeV2, 4 = 1.25 GeV^), B8 (m^ = 4.8 GeV, 
4 = 34GeV2, 4 = 1.5 GeV^). (a) ylf""''(0), (b) ^^"""(0), (c) ^^-^''(0). 

FIG. 7. The formfactors of D ^ p, normalized to their values at t = 0, as functions of t for all 
parameter sets and = 3GeV^. A pole- fit is sensible only for the vector formfactor and yields 
rupoi = (2.5 ±0.2) GeV. t^ax = 1.21 GeV^. (a) A^^'{t)/A^^'{0), (b) A^^''{t)/A^^'{0), (c) 

v^^p{t)/v^^p{o). 

FIG. 8. The formfactors oi B ^ p, normalized to their values at i = 0, as functions of t for 
all parameter sets and = 8GeV^. A pole-fit is sensible only for the vector formfactor and 
yields nipoi = (6.6 ± 0.6) GeV. tmax = 20.3 GeV^. (a) Af-^P{t)/Af^P{0), (b) A^-^f{t)/A^-'f{0), 
(c) V^-'P{t)/V^^P{0). 
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FIG. 9. The electron spectrum dT/dE for D — > pev as function of the electron energy E (set 
C5, M2 = 3GeV2). 

FIG. 10. (a) The electron spectrum dT/dE for B pev as function of the electron energy E 
(set B8, = 8GeV^). (b) Comparison of the spectra 1/T dT/dE as functions of E as obtained 
in this paper (parameters like (a)) and in the BWS |^ and the GISW model ||8|. The chosen 
normalization emphasizes the shape of the spectra. 

FIG. 11. Comparison of the spectra dT/dE as functions of E of the exclusive decays i? — > vr, pei^ 
(parameters like in Figs. ^ and |l^(a)) to the spectrum of the inclusive decay B — > X^eu, taken from 
pl| and calculated with the same parameters, (a) restsystem of the decaying B, (b) labsystem of 
a collider working at the T{4S) resonance. 

FIG. 12. Diagrams contributing to the Wilson-coefficient of the gluon condensate. Lines with 
a cross denote vacuum expectation values. 0, x, y are space-time coordinates, b, u, q denote quark 
flavours, q being a light quark (n or d). The weak vertex is at y. 
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